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The number of steps required in order to maximize a Bell inequality for arbitrary number of 
qubits is shown to grow exponentially with either the number of steps and the number of parties 
involved. The proof that the optimization of such correlation measure is a NP-problem is based 
on an operational perspective involving a Turing machine, which follows a general algorithm. The 
implications for the computability of the so called nonlocality for any number of qubits is similar to 
recent results involving entanglement or similar quantum correlation-based measures. 
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I. INTRODUCTION 

Quantum correlations lie at the heart of quantum in¬ 
formation theory. They are responsible for some tasks 
that posses no classical counterpart. It is plain from 
the fact that quantum measures are essential feature for 
quantum computation or secure quantum communica¬ 
tion, that one has to be able to develop some procedures 
(physical or purely mathematical in origin) so as to as¬ 
certain whether the state p representing the physical sys¬ 
tem under consideration is appropriate for developing a 
given non-classical task. Among those correlations, en¬ 
tanglement is perhaps one of the most fundamental and 
non-classical features exhibited by quantum systems PQ, 
that lies at the basis of some of the most important pro¬ 
cesses studied by quantum information theory m such 
as quantum cryptographic key distribution pj, quantum 
teleportation [7], superdense coding [8], and quantum 
computation mm- 

Other measures have been introduced in the literature 
that grasp features that are not captured by entangle¬ 
ment. They are not directly related to entanglement, 
but in some cases -specially when dealing with systems 
of qubits greater that two- they provide a satisfactory ap¬ 
proximate answer, like the maximum violation of a Bell 
inequality, that is, nonlocality. Local Variable Models 
(LVM) cannot exhibit arbitrary correlations. Mathemat¬ 
ically, the conditions these correlations must obey can 
always be written as inequalities -the Bell inequalities- 
satisfied for the joint probabilities of outcomes. We say 
that a quantum state p is nonlocal if and only if there 
are measurements on p that produce a correlation that 
violates a Bell inequality. 

Later work by Zurek and Ollivier m established that 
not even entanglement captures all aspects of quantum 
correlations. These authors introduced an information- 
theoretical measure, quantum discord, that corresponds 
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to a new facet of the “quantumness” that arises even for 
non-entangled states. Indeed, it turned out that the vast 
majority of quantum states exhibit a finite amount of 
quantum discord. Besides its intrinsic conceptual inter¬ 
est, the study of quantum discord may also have tech¬ 
nological implications: examples of improved quantum 
computing tasks that take advantage of quantum correla¬ 
tions but do not rely on entanglement have been reported 
[see for instance, among a quite extensive references-list 
P3HI5I]- Actually, in some cases entangled states are use¬ 
ful to solve a problem if and only if they violate a Bell 
inequality |16| . Moreover, there are important instances 
of non-classical information tasks that are based directly 
upon non-locality, with no explicit reference to the quan¬ 
tum mechanical formalism or to the associated concept 
of entanglement HZ]. A recent work studying how entan¬ 
glement can be estimated from a Bell inequality violation 
also sheds new light on the use of Bell inequalities [T5| 

In any case, the study of entanglement in multipartite 
quantum systems has been limited to few cases. As a 
consequence, other measures have been introduced in or¬ 
der the describe the “quantumness” of a certain (usually 
mixed) state p. In recent years the use of the maximum 
violation of a Bell inequality serves the purpose of de¬ 
scribing how nonlocal the state of the system is (see [19] 
and references therein). Although there is little connec¬ 
tion between entanglement and nonlocality (the former 
is based on how the tensor structure of the concomitant 
Hilbert space is split, whereas the latter ascertains how 
well a LVM can mimic quantum mechanics), nonlocality 
is a good candidate for describing correlations in quan¬ 
tum systems. 

To be more precise, the maximum violation of a Bell 
inequality for N parties B^ ax is the quantity chosen to 
approach entanglement is those scenarios [20] [21]. Thus, 
to know whether the computation of the maximum value 
of a Bell inequality is NP-hard seems a relevant and rea¬ 
sonable question. Previous approaches in the literature 
have dealt with the simplest possible instance of two par¬ 
ties [22], Alice and Bob, each possessing two nearly di- 
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chotomic observables. In the case of the Clauser-Horne- 
Shimony-Holt Bell inequality (CHSH) [23] , which is the 
strongest possible inequality for two parties (two qubits), 
it was proved that its maximum violation requires a com¬ 
putational work which grows exponentially with the num¬ 
ber of steps required, that is, it is a NP-problem [ 22 ]. 

In addition, a quantum measure such as discord has 
been recently proved to be NP-complete for the case of 
two qubits [2J. The fact that the computation of some 
entanglement and correlated quantities is NP-hard is usu¬ 
ally a consequence of the optimization involved in the 
definitions. The traditional tools required for optimizing 
Bell inequalities are borrowed from linear programming: 
inequalities are translated into convex polytopes, usually 
in high dimensions, and the proof for the general case 
of N parties involves a gigantic task which has not been 
successfully solved to date. 

The purpose of the present work is to provide an oper¬ 
ational approach to the process of carrying out the maxi¬ 
mization of a Bell inequality, based on a Turing machine, 
which will prove to be a NP-problem. A key ingredient 
will be the fact that Bell inequalities (not for probabil¬ 
ities) possess a recursive expression when they are gen¬ 
eralized to N qubits. In Section II we review previous 
results for CHSH for two qubits. In Section III we in¬ 
troduce the structure of the Bell inequalities employed, 
as well as the algorithm that the Turing machine will 
perform in this scenario. Finally, some conclusions are 
drawn in Section IV. 


II. PREVIOUS RESULTS 


The first approach to a Bell inequality (the CHSH in 
this case) for two qubits was carried out by Pitowsky 
[22'j . He carries out several extremely interesting investi¬ 
gations concerning the foundations of quantum mechan¬ 
ics. He also brings together high level characterizations, 
in geometrical language, of allowed classical and quan¬ 
tum correlation patterns. 

By using “classical correlation polytopes”, he provides 
significant insight into familiar the CHSH Bell-type in¬ 
equalities. Pitowsky provides an algorithm for finding 
the set of “generalized Bell inequalities” corresponding 
to any particular choice of the vectors defining the hy¬ 
perplanes of a convex polytope. Afterwards he proves it 
to be a not efficient one, that is, NP-complete. Thus, 
already for N = 2 qubits, the procedure of obtaining 
the maximum violation of the CHSH Bell inequality is 
inefficient. 




FIG. 1. (Color online) Plot of the typical performance of a 
simulated annealing maximizing the CHSH Bell inequality for 
a random state p of two qubits vs the number of Monte Carlo 
steps (MC). Every MC step contains 1000 random runs. Fig. 
(a) depicts the evolution of the maximum of CHSH for the 
state ||00)(00| + ||01)(01| + ||11)(11|, which is analytic (2/3). 
Fig. (b) depicts the same quantity for a random state p. See 
text for details. 


III. THE TURING MACHINE AND THE 
GENERALIZED OPTIMIZATION OF THE BELL 
INEQUALITY 


A. Bell inequalities 

Most of our knowledge on Bell inequalities and their 
quantum mechanical violation is based on the CHSH 
inequality [23 . With two dichotomic observables per 
party, it is the simplest [25] (up to local symmetries) 
nontrivial Bell inequality for the bipartite case with bi¬ 
nary inputs and outcomes. Let A\ and A 2 be two pos¬ 
sible measurements on A side whose outcomes are aj £ 
{—1,+1}, and similarly for the B side. Mathematically, 
it can be shown that, following LVM, \Bcf/s H W\ = 
laq&i + ai &2 + a 2 6 i — a 2 6 2 | < 2. Since ai(&i) and 
02 ( 62 ) cannot be measured simultaneously, instead one 
estimates after randomly chosen measurements the aver¬ 
age value S™, = Ex^hshW^W = E{A U B{) + 
E(Ai, B 2 )+E(A 2 , Hi)— E(A 2l -B 2 ), where E(-) represents 
the expectation value. Therefore the CHSH inequality 
reduces to 


\Bb V HSH\ < 2 - ( 1 ) 

Quantum mechanically, since we are dealing with 
qubits, these observables reduce to Aj(Bj) = aj(bj) ■ er, 
where aj(bj) are unit vectors in R 3 and a = (a x , a y , cr z ) 
are the usual Pauli matrices. Therefore the quantal pre¬ 
diction for 0 reduces to the expectation value of the 
operator B C hsh 


Ai 0 + Aj 0 B 2 + A 2 0 Bi — A 2 0 B 2 . (2) 

Tsirelson showed J2S] that CHSH inequality (|T]) is maxi¬ 
mally violated by a multiplicative factor \/2 (Tsirelson’s 
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bound) on the basis of quantum mechanics. In fact, it 
is true that \Tt(pabBchsh)\ < 2\/2 for all observables 
Ai, A 2 , Bi, B 2i and all states pab- Increasing the size 
of Hilbert spaces on either A and B sides would not give 
any advantage in the violation of the CHSH inequalities. 
In general, it is not known how to calculate the best such 
bound for an arbitrary Bell inequality, although several 
techniques have been developed m- 

A good witness of useful correlations is, in many cases, 
the violation of a Bell inequality by a quantum state. 
Although it is known that the violation of an TV-particle 
Bell-like inequality of some sort by an TV-particle entan¬ 
gled state is not enough, per se , to prove genuine multi¬ 
partite non-locality, it is the only approximation left in 
practice. 

The first Bell inequality for N = 3 qubits was pro¬ 
vided by Mermin 1281 . The Mermin inequality reads as 
Tr(pB M ermin ) < 2, where B Me rmin is the Mermin oper¬ 
ator 


BMermin — B aia2fl3 B ai b 2 b 3 Bfr ia2 B 3 Bb 1 b 2CL3 , (3) 

with B uvw = u-er(g)v-< 7 Cg>w-er with a = ( cr x , a y , a z ) 
being the usual Pauli matrices, and aj and bj unit vectors 
in R 3 . Notice that the Mermin inequality is maximally 
violated by Greenberger-Horne-Zeilinger (GHZ) states. 

In the case of N = 4 qubits, the first Bell inequal¬ 
ity was derived by Mermin, Ardehali, Belinskii and 
Klyslrko (MABK) (29]. The MABK inequality reads as 
Tv(pBmabk) < 4, where Bmabk is the MABK operator 

Bun — -E>ih2 ^ Hh21 — -B1211 — -B2111 — -B1122 — B1212 
~ B 2 h 2 — £>1221 — 02121 — B 2211 + B 2222 + £>2221 

+ -B2212 + -B2122 -®1222, 

(4) 

with B uvwx = u-(T®v-<r(g>w-crig)x-er with er = ( cr x , a yi a z ) 
being the usual Pauli matrices. We shall define 

Ha/abk = max Tt(pBmabk) (5) 

a n b j 

as a measure for the nonlocality content for a given state 
p of four qubits, aj and bj are unit vectors in R 3 . MABK 
inequalities are such that they constitute extensions of 
the CHSH inequalities with the requirement that gener¬ 
alized GHZ states maximally violate them. 

The optimization HOI HQ is taken over the two 
observers’ settings {aj, bj}, which are real unit vec¬ 
tors in R 3 . We choose them to be of the form 
(sin cos <pk , sin 9 sin (j>k, cos 9k). With this parameteri¬ 
zation, the problem consists in finding the supremum of 
Tt{pBmabk) over the (N = 4) {k = 1 • • ■ 4iV} angles. 

In the case of multiqubit systems, one must instead 
use a generalization of the CHSH inequality to N qubits. 


MABK inequalities are of such nature that they consti¬ 
tute extensions of older inequalities. To concoct an ex¬ 
tension to the multipartite case, we shall introduce a re¬ 
cursive relation [30] that will allow for more parties. This 
is easily done by considering the operator 

Bn +1 oc [(.Bi + B x ) ® .B/v + (Bi — B-jJ ® B N \, (6) 

with Bjv being the Bell operator for N parties and B\ = 
v • er, with a = ( a x , a y , a z ) and v a real unit vector. The 
prime on the operator denotes the same expression but 
with all vectors exchanged. The concomitant maximum 
value 

B“ ax = max Tr{pB N ) (7) 

a j ’ b j 

will serve as a measure for the non-locality content of a 
given state p of N qubits if aj and bj are unit vectors 
in R 3 . The non-locality measure ([7]) is maximized by 
generalized GHZ states, 2 t being the corresponding 
maximum value. 

However, there exist other measures m such as the 
Svetlichny inequalities [32] which serve the same purpose, 
having a similar structure extended to the TV-partite sce¬ 
nario [33] [31]. They have been used in the literature as 
entanglement-like indicators m- 

B. The Turing machine 

A Turing machine [35] has an infinite one-dimensional 
tape divided into cells. Traditionally we think of the tape 
as being horizontal with the cells arranged in a left-right 
orientation. The machine has a read-write head which is 
scanning a single cell on the tape. This read-write head 
can move left and right along the tape to scan successive 
cells. A table of transition rules will serve as the program 
for the machine. 

In modern terms, the tape serves as the memory of 
the machine, while the read-write head is the memory 
bus through which data is accessed (and updated) by 
the machine. One very important aspect is that we 
shall rely on the Turing-computability of the cost func¬ 
tion that maximizes the Bell inequality given a state pn- 
As known, there exists an entire class of these problems 
which is termed NP-complete (non-deterministic poly¬ 
nomial time complete) because the computational effort 
used to find an exact solution increases exponentially as 
the total number of degrees of freedom of the problem 
rise. 

As a consequence, approximated or heuristic meth¬ 
ods are required in practice for further analysis. The 
most successful statistical method to date is the stochas¬ 
tic model of simulated annealing introduced by Kirk¬ 
patrick, Gelatt, and Vecchi 36], that is, the Metropo¬ 
lis Monte Carlo algorithm with a fixed temperature T 
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FIG. 2. (Color online) Fig. (a) Describes how a Turing ma¬ 
chine computes the maximum violation of a Bell inequality 
for N = 4 qubits within the framework of simulated anneal¬ 
ing optimization. The number of terms to be evaluated by the 
machines grows exponentially as 4 N ~ 2 . Fig. (b) Emphasizes 
the fact that as the temperature drops, the available set of 
4N = 16 variables reduces until some precision on the cost 
function is reached. See text for details. 


at each state of the annealing schedule. There exist 
other methods which are not of statistical nature, such 
as downhill/amoeba or gradient methods [37], which in¬ 
volve finite differences when considering the correspond¬ 
ing function -a Bell inequality in our case- in terms of 
all real variables involved. 


each party randomly, keeping in the memory that some 
vectors will have the same position in the next move for 
some parties. The temperature is high at T 0 , which im¬ 
plies in Fig. 2 (b) that the domain of possible values for 
the variables Q is broadly spread. Keep in mind, how¬ 
ever, that every angle is reduced as follows: {6,, —► 6i 
mod Tr,ipi — t ipi mod 2-7t}. The machine then moves to 
the next term and performs similar operations accord¬ 
ingly. Finishing one cycle means visiting one after the 
other the entire 4 4 ^ 2 = 16 sites. After that, the machine 
has to compute Tr (BmPn)^ which is the cost function. 
Then, it starts the cycle anew with a different tempera¬ 
ture Ti (we can choose the temperature to decrease like 
T(s) = Toe -As , with s being the number of runs). As the 
temperature drops, the domain f l s shrinks as depicted in 
Fig. 2 (b). Thus, at every cycle the range of possible 
values for the 47V = 8 variables continuously decreases 
until we reach a desired precision, that is, the algorithm 
terminates when some stopping criterion is met. 

The basic algorithm is shown below. 


1 pjv initial multiqubit (N) state given 

2 T «— T 0 

3 repeat until stopping criterion is met 

4 repeat 4^” 2 times 

5 orientate N unit vectors 

6 move to the next term within the set 

7 endrepeat 

8 super operator B jv <— add all terms 

9 calculate Tr (BnPn) 

10 T <— T 0 e~ Xs 

11 endrepeat 

12 return B™ ax = Tr (Bnpn) 

( 8 ) 


C. Results 


If we do not want to specify a method in solving the 
optimization, we can rewrite the algorithm as: 


In either case -statistical or gradient-type method- we 
can program the Turing machine in the same way, be¬ 
cause after all it will undergo a Hamiltonian cycle chang¬ 
ing the value of several parameters of the total function to 
be optimized and at every step of the procedure. There¬ 
fore, we choose a simulated annealing approach to the 
program. 

Regarding Bell inequalities, one can choose the MABK 
or the Svetlichny inequalities to maximize for a given 
state p. In either case, owing to ([6]), the number of in¬ 
dividual terms grow exponentially as 4 N ~ 2 , N being the 
number of qubits. Let us then take the MABK inequali¬ 
ties. 

In order to illustrate how the Turing machine works, let 
us have the N = 4 case. The total number of independent 
variables are 4iV, but what makes the computation hard 
is the number of constraints that we have. The situation 
is depicted in Fig. 2 (a). The Turing machine reaches 
the first term in the tensor expansion of the MABK in¬ 
equality. It is free to move in space the unit vector of 


1 piv initial multiqubit (JV) state given 

2 repeat until stopping criterion is met 

3 repeat 4^ 2 times 

4 perform operations on N parties' observables 

5 register in memory 

6 move to the next term within the set 

7 endrepeat 

8 super operator Bn <— add all terms 

9 calculate Tr(BNP n) 

10 endrepeat 

11 return B^ ax = Tr (Bn Pn) 

(9) 


Every cycle contains at least 4 N ~ 2 visits, and the best 
computation of line 9 in the previous algorithm for the 
Turing machine is of 0(d 3 ), where d is the dimension of 
the square matrices (d = 2 N in our case) being multi¬ 
plied [38]. Therefore, we have an undefined number of 
times (at least two) x4 W-2 x (2 JV ) 3 number of steps re¬ 
quired to obtain U™ ax . In other words, at least we require 
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0(2 5N ) steps to solve the problem which, in view of the 
aforementioned result, clearly becomes NP with increas¬ 
ing number of parties. This is precisely the desired out¬ 
come: the computation of the maximum value of a Bell 
inequality requires a computational effort which grows 
exponentially with the number of parties N involved. 


IV. CONCLUSIONS 

Based on the iterative structure of the extension of 
Bell inequalities to the multiqubit case, we have shown 
that the maximization of the usual Bell inequalities em¬ 
ployed in the literature (except the ones for probabil¬ 
ities, as in [39]). an operation performed by a Turing 
machine, constitutes a NP-problem. This results some¬ 


how express the fact that, regarding nonlocality as a 
good resource for quantifying quantum correlations other 
than entanglement, the concomitant optimization be¬ 
comes non-tractable for high number of qubits. Further¬ 
more, even the fact of cheeking the plain violation of an 
inequality for a state pw, which implies Tt(Bnpn), is of 
0( 2 3Ar ), that is, it is limited in practice to a few number 
of qubits. 
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